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Abstract 

Context. Modem cosmology relies on the assumption of large-scale isotropy and homogeneity of the Universe. However, locally the 
Universe is inhomogeneous and anisotropic. This raises the question of how local measurements (at the ~ 10^ Mpc scale) can be used 
to determine the global cosmological parameters (defined at the ~ 10"* Mpc scale)? 

Aims. We connect the questions of cosmological backreaction, cosmic averaging and the estimation of cosmological parameters and 
show how they relate to the problem of cosmic variance. 

Methods. We used Buchert's averaging formalism and determined a set of locally averaged cosmological parameters in the context 
of the flat A cold dark matter model. We calculated their ensemble means (i.e. their global value) and variances (i.e. their cosmic 
variance). We applied our results to typical survey geometries and focused on the study of the effects of local fluctuations of the 
curvature parameter. 

Results. We show that in the context of standard cosmology at large scales (larger than the homogeneity scale and in the linear regime), 
the question of cosmological backreaction and averaging can be reformulated as the question of cosmic variance. The cosmic variance 
is found to be highest in the curvature parameter. We propose to use the observed variance of cosmological parameters to measure the 
growth factor. 

Conclusions. Cosmological backreaction and averaging are real effects that have been measured already for a long time, e.g. by the 
fluctuations of the matter density contrast averaged over spheres of a certain radius. Backreaction and averaging effects from scales in 
the linear regime, as considered in this work, are shown to be important for the precise measurement of cosmological parameters. 
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1. Introduction 

How do inhomogeneities in the matter distribution of the 
Universe affect our conception of its expansion history and our 
ability to measure cosmological parameters? Typically, these 
measurements rely on the averaging of a large number of in- 
dividual observations. In an idealised situation we can think of 
them as volume averages. To give an example, the power spec- 
trum, which is the Fourier-transformed two-point coiTelation 
function, may be seen as a volume average with weight e'*^*. 
Measurements of the properties of the large-scale structure rely 
on the observation of large volumes that have been pushed for- 
ward to ever higher redshifts in the last decade, from the two- 
degree Field survey (2dF, Colless et al. 2001) over the Sloan 
Digital Sky Survey (SDSS, Abazajian et al. 2009) to the cur- 
rent WiggleZ (Drinkwater et al., 2010) and Baryon Oscillation 
Spectroscopic Survey (BOSS, Eisenstein et al. 201 1). 

A theorem by Buchert states that the evolution of any 
volume-averaged comoving domain of an arbitrary iiTota- 
tional dust Universe may be described by the equations 
of a Friedmann-Lemaitre-Robertson- Walker (FLRW) model, 
but driven by effective sources that encode inhomogeneities 
(Buchert, 2000, 2001). The consequences of this have been 
studied extensively in perturbation theory (Kolb et al., 2005; 
Li & Schwarz, 2007, 2008; Brown et al., 2009a,b; Larena, 
2009; Clarkson et al., 2009; Buchert et al., 2011) and for non- 
perturbative models (Buchert et al., 2006; Marra et al., 2007; 
Rasanen, 2008; Kainulainen & Marra, 2009; Roy & Buchert, 
2010). Apart from the ongoing debate to what extent the global 
evolution is modified through backreaction effects from small- 



scale inhomogeneities (Rasanen, 2004; Ishibashi & Wald, 2006; 
Kolb et al., 2006; Buchert, 2008; Buchert et al., 2009; Wiegand 
& Buchert, 2010), Li & Schwarz (2008) showed that the mea- 
surement of cosmological parameters is limited by uncertainties 
concerning the relation between observable locally and unob- 
servable globally averaged quantities. 

In contrast to the well-studied cosmic variance of the cosmic 
microwave background, which is most relevant at the largest an- 
gular scales, the theoretical limitation on our ability to predict 
observations at low redshift aiises not only from the fact that we 
observe only one Universe, but also from the fact that we sample 
a finite domain (much smaller than the Hubble volume). Both 
limitations contribute to the cosmic variance. This is different 
from the sampling variance caused by shot noise, i.e. the limita- 
tion of the sampling of a particular domain because of the finite 
number of supernovae (SN) or galaxies observed. In the era of 
precision cosmology, the errors caused by cosmic variance may 
become a major component of the error budget. 

The purpose of this work is to demonstrate that the questions 
of cosmic averaging, cosmological backreaction, and the prob- 
lems of cosmic variance and parameter estimation are closely 
linked. We demonstrate that cosmic variance is actually one of 
the aspects of cosmological backreaction. 

There have already been many studies on the effect of the 
local dumpiness on our ability to measure cosmological param- 
eters. However, they focused mainly on the fluctuations in the 
matter density and on the variance of the Hubble rate. The for- 
mer question has gained renewed interest in view of deep red- 
shift surveys such as GOODS (Giavalisco et al., 2004), GEMS 
(Rix et al., 2004) or COSMOS (Scoville et al., 2007). Because 
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the considered survey fields are small, the variance of the matter 
density is an important ingredient in the error budget and it has 
been found to be in the range of 20% and more, demonstrated 
empirically in SDSS data by Driver & Robotham (2010) and 
calculated numerically from linear perturbation theory in Moster 
et al. (2010). 

The variance of the Hubble rate has been considered in the 
setup of this work, i.e. first order linear perturbation theory in 
comoving synchronous gauge, in Li & Schwarz (2008). A calcu- 
lation of the same effect in Newtonian gauge has been performed 
in Clarkson et al. (2009) and Umeh et al. (2010); the two meth- 
ods agree. Calculations of the fluctuations in the Hubble rate ow- 
ing to peculiar velocities alone have a longer history (Kaiser, 
1988; Turner et al., 1992; Shi et al., 1996; Wang et al., 1998). 
The effects of inhomogeneities on other cosmic parameters have 
been studied in Newtonian cosmology (Buchert & Ehlers, 1997; 
Buchert et al., 2000). 

There has been less activity in studying the effects of the 
third important player besides number density and Hubble ex- 
pansion: cosmic curvature. Even if the Universe is spatially flat, 
as expected by the scenario of cosmological inflation and con- 
sistent with the observations of the temperature anisotropics of 
the cosmic microwave background, the local curvature may be 
quite different. To answer the question how big this difference 
actually is for realistic survey volumes, we here extend the anal- 
ysis of Li & Schwarz (2007, 2008), where these effects have 
been estimated for the first time. For the spatially flat Einstein- 
de Sitter (EdS) model, the averaged curvature parameter Q.^ has 
been shown to deviate from zero by ~ 0.1 on domains at the 
100 Mpc scale. 

Here we adapt the analysis of Li & Schwarz (2008) to the 
case of a ACDM universe, introduce a more realistic power spec- 
trum and use observationally interesting window functions, not 
restricted to full-sky measurements. 

There has been quite some confusion about the choice of 
gauge and the dependence of the averaged quantities on it. 
Recently it has been shown (Gasperini et al., 2010) that this is 
not a problem if one consistently works in one gauge and then 
expresses the quantity that is finally observed in this frame as 
well. This is easier in some gauges than in others, but the result 
is (as expected) the same, as is also confirmed explicitly by the 
fact that our results are consistent with those of Clarkson et al. 
(2009) and Umeh et al. (2010), obtained in Newtonian gauge. 
The quest for simplicity explains our choice to use comoving 
synchronous gauge, because this is the frame that is closest to 
the one used by the observers. 

Fig. 1 depicts the theorist's and the observer's view on the 
Universe in a schematical way. It points out that in the end it 
is the average quantities that we are interested in, but that in an 
intermediate step, observers like to think of the objects they mea- 
sure to lie in a comoving space with simple Euclidian distances. 
The comoving synchronous gauge, which has a clear notion of 
"today" for a fluid observer sitting in a galaxy as we do, helps to 
define the things observers measure in a simple way. 

In Section 2 we establish the conceptual framework for the 
study of effects of inhomogeneities on observable quantities. 
Section 3 and 4 generalize some of the results of Li & Schwarz 
(2007, 2008) and Li (2008) from an Einstein-de Sitter (EdS) 
to a ACDM background and implement a more realistic mat- 
ter power spectrum. Section 5 investigates the effects of various 
window functions, again extending the analysis of Li & Schwarz 
(2008); Li (2008). Section 6 concentrates on deriving the mag- 
nitude of curvature fluctuations for realistic window functions. 
Section 7 applies the formalism to the local distance measure 
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Figure 1. Comparison of the theorist's and observer's view on 
the Universe. Our calculation in comoving synchronous gauge 
facilitates the description of the boundaries of the experimen- 
tally investigated regions in our Universe. Note, that recently 
there have been attempts by Bonvin & Durrer (201 1) to directly 
relate the two upper circles. This was done by calculating the 
predictions for the quantities in redshift space explicitly from 
the perturbed ACDM model. 

Dy, determined in the observation of baryonic acoustic oscil- 
lations (BAO). Section 8 is a remark on the link between the 
variance of averaged expansion rates at different epochs and the 
background evolution, before we conclude in Section 9. 

2. Inhomogeneity and expansion 

We assume that the overall evolution of the Universe is described 
by a flat ACDM model, which we adopt as our background 
model throughout this work. Global spatial flatness does not pre- 
vent the local curvature to deviate from zero. 

The distribution of nearby galaxies indicates that the scales 
at which the Universe is inhomogeneous reach out to at least 
100 Mpc. Above these scales it is not yet established if there 
is a turnover to homogeneity, as was claimed in Hogg et al. 
(2005), or if the correlations in the matter distribution merely 
become weaker, but persist up to larger scales, as is discussed 
in Labini (2010). At least morphologically, homogeneity has not 
been found up to scales of about 200 Mpc (Kerscher et al., 1998, 
2001;Hikage et al., 2003). 
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Consequently, we need a formalism that is applicable in the 
presence of inhomogeneities at least for the description of the 
local expansion. This may be accomplished considering spa- 
tial domains D and averaging over their locally inhomogeneous 
observables (Buchert, 2000, 2001). Technically one performs 
a 3+1 split of spacetime. Because we considered pressureless 
matter only, we chose a comoving foliation in which the spa- 
tial hypersurfaces are orthogonal to the cosmic time. This means 
that the formalism will not be able to take into account light- 
cone effects. Therefore, it is well adapted for regions of the uni- 
verse that are not too extended, a notion described more pre- 
cisely in Sect. 5. The equations then describe the evolution of 
the volume of the domain £), given by |£)|^ := dyUg, where 

djUy := [^^^^(f, x)]'^^d-';ic and ^^^g is the fully inhomogeneous 
three-metric of a spatial slice. 

To obtain an analogy to the standard Friedmann equations, 
one defines an average scale factor from this volume 



(1) 



where the subindex denotes "today", as throughout the rest 
of this work. The definition implies = 1. In analogy to the 
background model i/© :- a^olciD- The foliation may be used to 
define an average over a three-scalar observable (9 on a domain 
D in the spatial hypersurface 



(0)»(f): = 



^ (9(f,x)d//g 



(2) 



Examples for these observables include the matter density g or 
the redshift z of a group of galaxies in the domain D. 

The scalar parts of Einstein's equations for an inhomoge- 
neous matter source become evolution equations for the volume 
scale factor driven by quantities determined by this average: 



1 1 



3/?^ = 8;rG <e>» - - <!R>» - -Q© + A 



(3) 

(4) 
(5) 



The expansion of the domain D is determined by the average 
matter density, the cosmological constant, the average intrinsic 
scalar curvature {H)^ and the kinematical backreaction Qj^. The 
latter encodes the departure of the domain from a homogeneous 
distribution and is a linear combination of the variance of the 
expansion rate and the variance of the shear scalar. 

Equations (3) to (5) mean that the local evolution of any in- 
homogeneous domain is described by a set of equations that cor- 
responds to the Friedmann equations. 

The cosmic parameters, defined by 



D ._ 
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(6) 
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are domain-dependent. Owing to the fluctuating matter density, 
the curvature and the average local expansion rate will also fluc- 
tuate. When we constrain ourselves to the perturbative regime, 
the modification due to Qj^ is important on scales of the order of 
10 Mpc (Li & Schwarz, 2008). 



Here we are interested to see to what extent the values of the 
parameters (6) vary if we look at different domains of size D in 
the Universe. We could therefore define an average over many 
equivalent domains D at different locations in our current spa- 
tial slice. For an ergodic process, however, this is the same as the 
variance of an ensemble average over many realizations of the 
Universe keeping the domain D fixed, but changing the initial 
conditions of the matter distribution. This is the quantity that we 
calculate in theory and therefore we have to rely on the assump- 
tion of ergodicity when comparing our results with the observa- 
tion. In our case this ensemble average is taken over quantities 
that are volume averages. This means that for any observable O 
there are two different averages involved. The domain averag- 
ing, (O)^, and the ensemble average, O. We assume that both 
averaging procedures commute. 

The fluctuations are then characterized by the variance with 
respect to the ensemble averaging process. 



'ri{0)o):^({0)l-{0}oy 



(7) 



An example for a common observable calculated in this man- 
ner would be cTg as the ensemble rm.s. fluctuation of the matter 
density field. To quantify these fluctuations of a general observ- 
able O, we use the theory of cosmological perturbations. In Li 
& Schwarz (2007) it has been shown that to linear order, the 
(3o-term in equations (4) and (3) vanishes, {g}^, {K}^ and //© 
however, have linear corrections. 

Furthermore, Li (2008) argued that there is no second-order 
contribution to the fluctuations for Gaussian density perturba- 
tions if their linear contributions are finite. This may be seen 
by decomposing the observable O into successive orders O = 
^(0) + ^(1) + Q{2) + ... It is usually assumed that = 0. Now, 
for Gaussian perturbations only terms of even order give rise to 
non-trivial contributions. Therefore, (7) may be expressed as 

x2 



o-(0) = ^fiO'^^) 



1 -t- 



2(0(i)r 



(8) 



which shows that the correction to the leading order linear term 
is already of third order This is why we content ourselves for the 
evaluation of the fluctuations in the parameters (p)^,, {'R}^) and 
Hj) or Q®, and //^ to a first order treatment. This argument 

does not apply if (OO)^ - 0, as is the case for (3©. In that case 
(3o and o-{Qo) are of second order in perturbation theory. 

3. Cosmological parameters and their mean from 
local averaging 

The analysis of this work is based on standard perturbation the- 
ory in comoving (synchronous) gauge and we use results and 
notation of Li & Schwarz (2007). The perturbed line element 

ds^ = fl2 |_^^2 ^ _ 2^(1)) ^. . + D..^Wj dx'dx>} (9) 

defines the metric potentials 1/^0(77, x) and;(f(')(77, x). Below we 
use the convention ao = 1 for today's scale factor. We use 
conformal time 77 and the traceless differential operator Djj = 
didj - \6ijA on a spatially flat background. The geometrical 
quantities of interest are the local expansion rate and the spa- 
tial curvature. The former follows from the expansion tensor and 
reads 



(I)' 



(10) 
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where ()' stands for the derivative with respect to conformal 
time. Calculating the spatial Ricci curvature from the above met- 
ric yields 



(11) 



By the covariant conservation of the energy momentum tensor, 
i/f*'* is related to the matter density contrast 



6(ri,x) :-- 



P 



9(0) 



by 



1 



(12) 



(13) 



demonstrating that the importance of curvature effects grows 
proportional to the formation of structures. A remarkable prop- 
erty is that 2 ^iifl,^) - 1 holds not only for the FLRW back- 
ground, but also at the level of perturbations. For linear pertur- 
bations the kinematic backreaction term does not play any role, 
but becomes important as soon as quadratic terms are consid- 
ered. 

Let us now compare these local quantities with the domain- 
averaged expansion rate and the spatial curvature (Li & Schwarz, 
2008). From the definition of the average ()^ in (2) we find that, 
in principle, fluctuations in the volume element d^g have to be 
taken into account. Writing dju^ — Jd^x with the functional de- 
terminant J - {\ - 3i/'*''), the average over the perturbed hy- 
persurface agrees with an average over an unperturbed Euclidean 
domain 



with I (x) denoting a constant of integration. This constant plays 
no role in the following, because 6 and "R involve only time 
derivatives of 

For dust and a cosmological constant, Einstein's equations 
give the well-known relation 



/oO(»7dx /gO'Ddx 



(24) 



„ a' , ^nGpf 

5" + -5' = —6, 

a a 



(14) 



if we restrict our attention to linear perturbations. 

We express domain-averaged quantities in terms of the vol- 
ume scale factor aj), because we assume that the measured red- 
shift in an inhomogeneous universe is related to the average scale 
factor. This has been advocated by Rasanen (2009), where the 
relation 



at first order in perturbation theory. For the ACDM model the 
solution reads [see, e.g. Green et al. (2005)] 



(1+z) 



(25) 



= z)(T) 



(15) 



where 6q (x) is the density perturbation today. D (a) is the growth 
factor given by 



111 

D(a) = 02/^1 1 1, "^'"'^'^ 



with c = 



(16) 



has been established. Note that in principle one would have to 
introduce averaging on some larger scale than D to connect this 
background average on some domain S to the redshift. For the 
sake of simplicity, and because we content ourselves with small 
redshifts, we use the same domain D. This limits the validity of 
the result to small redshifts. 

In order to relate a and a©, we start from 



1 fljj I I a' 



and 2F1 is a hypergeometric function. In the following we denote 
today's value of the growth factor by Do = D(l). 

Plugging this solution into (10) and using (13), we find the 
local expansion rate 



1 dn 

3 a-D a a-D a \ a 

To first order this relation gives 
1 



all 



(26) 



(27) 



^e(a, x)^Ho^J^ VTT^ |l - ^/ (a) 6 (a, x)j , (17) We finally obtain the averaged Hubble rate 



expressed in terms of the growth rate 
dlnD(a) 



f{a):-- 



5— — 3 

D(a) 




1 - 



D(an) 



-3 



Do 



6(1 + cfl^) 



Da 



(28) 



dlno 2(1-1- ca^) 
From (1 1) we find the local spatial curvature 



' Do ■ 



From these quantities we can define local Q. functions. 



1 



1 + Cfl3 

1 



(a,x) 
(a, x) 
^^A(fl,x) = 

(1 + ca^) 
QQ(fl,x) = 0, 



l+|l + -/(fl)|5(fl,x) 



1 + ca 



3 + 3/(0) 



6 {a, x) , 



^3 



1 + -/(fl)5(fl,x) 



(18) 

(19) 

(20) 

(21) 

(22) 
(23) 



and the averaged spatial curvature 

^0 (^0)£> 



00. 



10 Q„ 



aj, Dq 



For later convenience we also define the function 



fviav) : = 



5^^ - 3 

Do 

2(1 -i-cfl^)' 



(29) 



(30) 



which is our modified version of the growth rate of Eq. (18), 
multiplied by D (a) /Dq. It basically encodes the deviation of the 
time evolution of the Hubble perturbation from the time evolu- 
tion of the ensemble averaged Hubble rate 



Ho (an) 




(31) 
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as may be seen from the resulting expression 
1 



Ho = Ho (a©) 1 - -fo (fl©) <5()>a 



In the Einstein-de Sitter limit (c — > and Q.„ 



at 



Ho = (l - ^fl0 {5())v 



(32) 
1) we arrive 

(33) 
(34) 



In order to compare this with the results of Li & Schwarz (2007, 
2008), we define the peculiar gravitational potential ip (x) via 

2 a 3 tia 



(35) 



and obtain 



Ho - 
and 

Wo 



3fo 



-3/2 

7 



1 



20 



-2 



(36) 



(37) 



While our results agree for the spatial curvature, is differ- 
ent from the result in Li & Schwarz (2008), because there the 
assumption a I was made when applying (27). 

Let us now turn to the dimensionless Q®-parameters. To first 
order, they may be expressed as 



1 



1 + ca^ 



1+11 + 3 ./o (flo) I <^())£> 



^ + (a») 

I + cat, J 



(<5o)© , 



1 + ca^ 



1 + 



(39) 



(40) 



(41) 



When taking the limit 2) — > in Eqs. (38) - (41), we recover 
the point-wise defined O-parameters of Eqs. (20) - (23). This 
provides a self-consistency check of the averaging framework. 

From the expressions for the Q®-parameters one can eas- 
ily calculate the ensemble averages and the ensemble variance. 
{^o)d - 0' since the domain-averaged overdensity of £), in gen- 
eral non-zero, averages out when we consider a large number of 
domains of given size and local density fluctuations drawn from 
the same (Gaussian) distribution. 

Here we adopt the common view that linear theory is a 
good description of the present universe at the largest observable 
scales (which has been questioned recently in Rasanen 2010). 
We then find the ensemble average of the curvature parameter 

to vanish. For the matter density parameter Eq. (38) yields 



^Siao) = (l +cfl^) 



(42) 



This may be used to verify that the relation + - 1 holds. 

In addition, this relation implies that («£>„) corresponds to 
today's background matter density parameter: 



(43) 



However, this is true at first order in the density contrast only, 
because in this case ensemble averages agree with background 
quantities. At higher orders, the ensemble averages diff'er from 
the background quantities. 



4. Variances of locally averaged cosmological 
parameters 

After having convinced ourself that the expectations of the av- 
eraged Q®-parameters are identical to their ACDM background 
values up to second-order corrections, we now turn to the study 
of their ensemble variances. 

All variances of domain averaged cosmological parameters 
can be related to the variance of the overdensity of the matter 
distribution, o-({6o)o)- 

In order to specify (5o)o, we introduce the normalized win- 
dow function (X) and write 



0)0 - ^0 

-[ 

Jk3 



^o(k)W»(k)d3fe, 



(44) 



A tilde denotes a Fourier-transformed quantity. With the defini- 
tion of the matter power spectrum 



5o(k)5o(k') = 5""''(k + k') Poik), 



(45) 



where S^""^^ denotes Dirac's delta function, the ensemble vari- 
ance of the matter overdensity becomes 



(38) {crof -.^ cr^ {{SoW 



(/t)W2)(k)W2)(-k)d^L (46) 



For a spherical window function, this expression is the well- 
known matter variation in a sphere, often used to normalize the 
matter power spectrum by fixing its value for a sphere with a 
radius of 8/;"'Mpc (cr^). To calculate this variance, we assume 
a standard ACDM power spectrum in the parametrization of 
Eisenstein & Hu (1998). Knowing cro at a particular epoch of 
interest, we can calculate all the fluctuations in the cosmic pa- 
rameters. They read: 



1- 



6Hd = -//©(fl©)/©(fl©)cro^ , 



6nl 



3 

of (flo) 



1 + ^/. 



1 4- 



2 fv (ao) 



(47) 
(48) 

(49) 



6^2 



(50) 

0{icrof) , (51) 

where e.g. (5f2® denotes the square root of the variance, : = 
cr^Q®^ Q®, and (a^) were defined in Eq. (42), (31) and 

(30)andQf = l-Of. 

These variances are the minimum ones that one can hope to 
obtain by measurements of regions of the universe of size D. 
They do not include any observational uncertainties, nor biasing 
or sampling issues. They are intrinsic to the inhomogeneous dark 
matter distribution that governs the evolution of the Universe. 



5 



A. Wiegand, D.J. Schwarz: Inhomogeneity-induced variance of cosmological parameters 




Figure 2. The two survey geometries considered (separately). A 
simple cone with one single opening angle a and a slice given 
by two angles p and y. 



Equations (47) to (51) are interesting in two respects: Firstly, 
our expression for 6Hj] is simpler than the one in Umeh et al. 
(2010), nevertheless, both results agree. Secondly, Eqs. (47) to 
(51) quantify the connection between fluctuations in cosmologi- 
cal parameters and inhomogeneities in the distribution of matter 
If we choose "today" as our reference value, Eqs. (47) to (50) al- 
low us to predict the domain averaged cosmological parameters: 



Ho - Hq 
Q| = 

with 



+ |/0o) 



(52) 



/B„=/B(fl©„)=^(5Z)o'-3) 



0.5 ACDM 
1.0 EdS 



(53) 



where we assumed f2,„ = 0.3 for ACDM. More generally, for 
f2„, > 0.1, /©„ may be approximated by (Lahav et al., 1991; 
Eisenstein & Hu, 1998) 



(54) 



From the knowledge of cr^^ we may therefore easily derive 
the variation of cosmological parameters. To relate our calcula- 
tions to real surveys, we elaborate in the next section on how to 
calculate ctj;)^ for several survey geometries. 

5. The effect of the survey geometry 

Observations of the universe are rarely full-sky measurements 
and typically sample domains much smaller than the Hubble 
volume. We therefore must address the problem of the survey 
geometry. Effects from a limited survey size are in particular im- 
portant for deep fields, as studied for example in Moster et al. 
(2010). While in their case, for small angles and deep surveys, 
approximating the observed volume by a rectangular geometry is 
appropriate, it probably is not appropriate for the bigger survey 
volumes that we have in mind. 

We therefore chose two different geometries that resemble 
observationally relevant ones. Firstly, we used a simple cone 
with a single opening angle a. The second geometry is a slice 
described by two angles [5 and y for the size in right ascension 




Volume (/tq yMpc') 



Figures. Variance of the matter density, cr^,, as a function of 
the observed domain volume. Data were derived from the SDSS 
main sample by Driver & Robotham (2010). The dashed (red) 
line shows the fit of Driver &. Robotham (2010) to the data, 
the solid (green) line is our result including the sample variance 
[solid (blue) line at the bottom]. 



and declination respectively. In the radial direction we assumed 
a top hat window, whose cut-ofF value corresponds to the depth 
of the survey. Both shapes are shown in Fig. 2. 

To calculate cr^Oa for both geometries, we used a decompo- 
sition into spherical harmonics. This allowed us to derive an ex- 
pression for the expansion coefficients in terms of a series in 
cos {2na) for the cone and a similar one for the slice, depending 
on trigonometric functions of [5 and y. The radial coefficients 
were calculated numerically using the ACDM power spectrum 
of Eisenstein & Hu (1998), including the effect of baryons on 
the overall shape and amplitude of the matter power spectrum, 
but without baryon acoustic oscillations. 

All plots use best-fit ACDM values as given in Komatsu et al. 
(2011); Oi = 0.0456, Qedm = 0.227 and = 0.963. The power 
spectrum is normalized to erg = 0.809. 

To ensure that the result of our calculation for the slice- 
like geometry and a standard ACDM power spectrum is reason- 
able, we compared it with an analysis of SDSS data by Driver 
& Robotham (2010). In Fig. 3 we show this comparison of 
their rm.s. matter overdensity cr^, obtained from the SDSS main 
galaxy sample, in terms of its angular extension (and hence the 
volume). The dashed line going through the points shows their 
empirical fit to the data. The solid green line shows our result for 
the cosmic variance of a sUce with respective angular extension 
(for y6 = y), plus their sample variance. Note that our result is 
not a fit to SDSS data, but is a prediction based on the WMAP 
7yr data analysis. Additionally, the real SDSS window function 
is slightly more complicated than our simplistic window, thus 
perfect agreement is not to be expected. 

For the full SDSS volume, cr^,^ is shown in Fig. 4. For com- 
parison we also added the smaller, southern hemisphere 2dF 
survey and a hypothetical full sky survey. For the two surveys, 
we assumed an approximate angular extension of 120°x60° for 
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Figure 4. Variance of the matter density, cr^), for survey geome- 
tries resembling the 2dFGRS, the SDSS and a hypothetical full 
sky survey as a function of maximum redshift considered. We 
find that the determination of the local crj] below redshifts of 
0. 1 (corresponding to ~ 400 Mpc) is fundamentally limited by 
cosmic variance to the 1% level. 



SDSS and for the two fields of the 2dF survey 80°xl5° and 
75°xl0°. The ongoing BOSS survey corresponds to the plot for 
the SDSS geometry because it will basically have the same an- 
gular extension. Because it will target higher redshifts, it is not 
in the range of our calculation, however. As a rough statement 
(the precise value depends on the redshift), one may say that 
the 2dF survey is a factor of 5 and the SDSS survey a factor of 
2.5 above the variance of a full sky survey. This is interesting 
because the SDSS survey covers approximately only 1/6 of the 
full sky and the 2dF survey only 1 /20. This is due to the angular 
dependence of crj). We find that fluctuations drop quickly as we 
increase small angles and flattens at large angles. 

From Fig. (4) we see that the cosmic variance of the matter 
density for the SDSS geometry is 5% at a depth of z ^ 0.08 and 
still 1% out to z ~ 0.23. Note, however, that the extension of 
the domain of (spatial) averaging to a redshift of 0.35 is clearly 
not very realistic because lightcone effects will become relevant 
with increasing extension of the domain. The assumption that 
this domain would be representative for a part of the hypersur- 
face of constant cosmic time becomes questionable. We expect, 
however, that in this range evolution effects will only be a minor 
correction to the result presented here because of the following 
estimation: 

The main effect we miss by approximating the lightcone by 
a fixed spatial hypersurface is evolution in the matter density 
distribution. To specify what "region that is not too extended 
spatially" means, one should therefore estimate the maximum 
evolution in a given sample. This can be done by determining 
the growth of the density contrast in the outermost (and there- 
fore oldest) regions of the sample. In the linear regime con- 
sidered here, the evolution of the density contrast is given by 
6(z) - 6oD{0) /D{z). Therefore, lightcone corrections should 
be smaller than e, = 1 - D (z) /D (0) which is e, « 5% for z = 0. 1 
and ei * 14% for z = 0.3. Therefore, the order of magnitude 
of our results should be correct on a wider range of scales, but 
on scales above z = 0.3 the corrections to our calculation are 
expected to pass beyond 15%. 

Finally, it should be noted that for large volumes the actual 
shape of the survey geometry is not very important. As long 
as all dimensions are bigger than the scale of the turnover of 



the power spectrum, the deviation of the cosmic variance for 
our shapes, compared to those of a box of equal volume, is at 
the percent level. To reach this result, we compared crj]^ for 
the slice-like geometry to its value for a rectangular box of the 
same volume. We used a slice for which fi - j. The box was 
constructed to have a quadratic basis and the same depth as the 
slice in radial direction. Therefore the base square of the box 
is smaller than the square given by the two angles of the slice. 
The result of this comparison is that the deviation of cr^Jj' from 
the value for the slice is at most 6% for angles above /? ^ 10°. 
For smaller angles the deviation becomes bigger and redshift- 
dependent. This is caused by the changing shape of the power 
spectrum at small scales. The large angle behavior confirms an 
observation of Driver & Robotham (2010). They found that the 
cosmic variance in the SDSS dataset was the same for both of 
the two geometries they considered. 

6. Fluctuations of the curvature parameter 

After the general study of the effect of the shape of the observa- 
tional domain D on ctj)^, one may ask for which parameter the 
fluctuations are most important. 

The three lowest lines in the plot of Fig. 5 show that this 
is the case for the curvature fluctuations. The two lowest lines, 
showing the fluctuations 6Q.„, and ^Z/©,, IH^Oo for the full sphere, 
lie a factor of 1.6 and 3.8, respectively, below the respective cur- 
vature fluctuations SQ.^. Therefore the fluctuations of Q„, play 
a smaller role for all universes with Q„, < 1 . The uncertainty in 
//o,,, which has been in the focus of the investigations so far (Shi 
et al., 1996; Li & Schwarz, 2008; Umeh et al., 2010), contributes 
even less to the distortion of the geometry, as we shall discuss in 
Section 7. 

What this means for real surveys, such as the 2dF or the 
SDSS survey, is shown by the three upper lines in Fig. 5. They 
compare 6Q.^ for the slices observed by these surveys to that 
of a full sky measurement. 6Q.^ is bigger than one percent up 
to a redshift of 0.18 for the SDSS and 0.28 for the 2dF survey 
and it does not drop below 0.001 for values of z as high as 0.5. 
This may seem very low, but it has been shown that getting the 
curvature of the universe wrong by l%c already affects our abil- 
ity to measure the dark energy equation of state w (z) (Clarkson 
et al., 2007). Of course one has to keep in mind that for high 
redshifts one has to be careful with the values presented here be- 
cause they are based on the assumption that the observed region 
lies on one single spatial hypersurface. Because this approxima- 
tion worsenes beyond a redshift of 0.1, there may be additional 
corrections to the size of the fluctuations stemming from light- 
cone effects. 

To investigate the curvature fluctuations for more general ge- 
ometries, we show in Fig. 6 the angular and radial dependence 
of the curvature fluctuation 5D.^ (a) for the cone-like window of 
Fig. 2. 

On the l.h.s. of Fig. 6 we evaluate the angular dependence. 
For a survey that only reaches a redshift of 0.1, the fluctuations 
are still higher than 0.01 for a half-sky survey. It is interesting to 
note that for a deeper survey, 5Q.^ {a) grows much faster when 
a is reduced than for a shallow survey. This is because cr^^ (R) 
changes from a relatively weak 7?"' decay to a R^^ decay on 
larger scales. For z = 0.35, this behavior dominates and a de- 
crease in a increases cr^^ (R, a) stronger than in the R^^ regime. 

On the r.h.s. of Fig. 6 we show the dependence of 6Q.^ 
on the survey depth for some opening angles of the cone-like 
window. For narrow windows the fluctuation in stays high. 
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Figure 6. Cosmic variance of the curvature parameter. Left panel: Dependence of on the opening angle of the cone-like survey 
geometry of Fig. 2 for different top hat depths of the survey. Right panel: Dependence of 6Q.^ on the depth of the survey. For a 
small cone of 6° opening angle we expect curvature fluctuations of 10% up to 450/z 'Mpc. 
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Figures. Top three lines: the expected r.m.s. fluctuation of the 
curvature parameter, for geometries resembling the 2dF- 

GRS, the SDSS, and a full sky. The two lowest lines are the 
expected r.m.s. fluctuations of the parameters Q^" and Hj)^ for a 
full-sky survey extending to the respective redshift. The curva- 
ture fluctuations turn out to be higher than all other fluctuations. 



even beyond the expected homogeneity scale of lOO/i 'Mpc. For 
R - 200/i"'Mpc and a 6° window, for example, it is still at 
6Q.^ ~ 0.2. For smaller beams these fluctuations persist even out 
to much longer distances. Therefore they play an important role 
for deep field galaxy surveys, as shown in Moster et al. (2010); 
Driver & Robotham (2010) for the matter density fluctuations. 
But even for wider angles, fluctuations in curvature persist on 
sizeable domains. If one recalls that the distance given for the 
full sphere of 360° is its radius, this means that regions in the 
Universe as big as 540/z^'Mpc have typical curvature fluctua- 
tions on the order of 1%. This is not that small because the last 
scattering surface ai z ~ 1100 is only 9600/;"'Mpc away. One 
of these regions therefore fills more than 5% of the way to that 
surface. 

To put these values into perspective, we compare the 
WMAP Syr confidence contours (Komatsu et al., 2009) on the 
curvature parameter with those that may in principle be derived 
from the 2dF or the SDSS survey in Fig. 7. Because they only 



sample a finite size of the Universe, one cannot be sure that this 
value is indeed the background value and not only a local fluc- 
tuation. The cosmic variance induced by this finite size effect 
is, for the 2dF survey volume up to z s; 0.2, shown by the two 
second-largest (red) ellipses. The two innermost (blue) ones de- 
pict the minimum possible error using the SDSS survey volume 
up to z « 0.3. Clearly, the determination of may perhaps be 
improved by a factor of two if one were to eliminate all other 
sources of uncertainty. This may be less if lightcone effects play 
a non-negligible role already for z ~ 0.3. 

Fig. 8 shows the dependence of the curvature fluctuation on 
the considered cosmology. For this study we fixed the spectral 
index and the normalization of the spectrum to = 0.963 and 
erg - 0.809, respectively. We varied each of the other parame- 
ters one after another, while keeping the remaining ones fixed 
to the concordance values. We used the SDSS geometry out to 
a redshift of z = 0.09 as a reference value at which we con- 
ducted this investigation, because the concordance values lead 
to a of 0.01 for this configuration. Interestingly enough, the 
dependence on the Q„, parameter is very weak. This means that 
the value does not differ much for the flat ACDM model and 
the EdS model. This is surprising, because the prefactor of cr^^ 
in (49) changes by a factor of 3 from about 5/1 1 for ACDM to 
5/3 for EdS. This rise, however, is compensated for by a drop of 
the value of cr^yg. The reason for this drop is that a higher Q,„ 
leads to more power on small scales. Because we kept the inte- 
grated normalization fixed at a given value of erg, this means less 
power on large scales, i.e. at z = 0.09. Moreover, a variation of 
the Hubble constant h and the baryon fraction fi, has only a small 
effect around the concordance value. changes significantly 
only for more extreme values of and h. 



7. Fluctuations of the acoustic scale 

Let us now turn to the effect of fluctuations caused by inhomo- 
geneities on the local distance estimates. An important distance 
measure, recently used in BAO experiments, is Dy- It was intro- 
duced in Eisenstein et al. (2005) and mixes the angular diameter 
distance and the comoving coordinate distance to the BAO ring. 
It is measured through the BAO radius perpendicular to the line 
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Figure 7. Minimum confidence contours in and achiev- 
able in different volumes through fluctuations of matter. The 
green (outermost) ellipses are the 95% and 68% contours for 
the volume from which the HST data are drawn. The next inner 
(red) ones are for a survey of the size of the 2dF survey up to 
z - 0.2. In the middle there is a small double ellipse in blue, 
showing the values for the SDSS volume up to z = 0.3. The 
background image depicts the results from WMAP 5 (Komatsu 
et al., 2009). They give the experimental values and uncertainties 
on these parameters for a combination of various experimental 
probes. 

of sight and the comoving radius parallel to the line of sight 
r\\. 



rba. 



(55) 



One can, therefore, determine the distance Dy to the correspond- 
ing redshift, if the comoving radius of the baryon ring ri,„o is 
known. This may be achieved by a measurement of the angle of 
the BAO ring on the sky AO and its longitudinal extension Az/z- 
The precise definition of Dy is derived from the expressions of 
the comoving distances ry and r^: 



z+Az 



ni 



J H(z') 



cAz 



cz Az 



H(z) H(z) z 



ri. = (1 +z)Da(z)A0, 
from which we find 



Dviz) = 



cz 
Hiz) 



Dliz) 



where Dm is the comoving angular distance 

-1 



Dm(z) = c ( ^/n,Hoy sinh ( (z)) , 
with 
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Figures. Dependence of 5f2® on some cosmological parame- 
ters for a spherical domain extending to z = 0.09. The basis is 
the ACDM model with = 0.0456, Q-cdm = 0.227, h = 0.7, 
«s - 0.963 and erg = 0.809. For this model and for the cho- 
sen redshift, ~ 0.01. We then varied Q,„ = -H ^cdm^ 
fb - Q.blO.„, and h between and 1, holding the other param- 
eters fixed at their aforementioned values. Because erg is fixed, 
the fluctuation in is nearly independent on 0,„. 



As already mentioned above, the term Q.^ vanishes in our 
first-order treatment and the curvature contribution scales as a^. 
Therefore we may express the Hubble rate as 

Hpjz) 



= [(1 + z)3q^° + (1 + z)2q^° + (1 - Q,^° - Q^")]^ , (61) 



where we assumed the relation between redshift and average 
scale factor of Eq. (25). We may now calculate the fluctuation 
ofDy, 



5r^ 



5r^ 
ri. 



dDv 



l-{l+zf Hi 



Hoizf 



l-(l+z)^ Hk 



5D^ 



2 



, l'(z) 



HoizY 



I'iz) 



KzY , 

6 Kz) 



6D:, 



Da 



Kz) 
1 ^ ^ 
3 ry 



377 



(62) 

(63) 
(64) 



where I'iz) denotes a partial derivation with respect to the re- 
spective parameter, i.e. Q^*" or 0^°. Note that I'iz) and Kz) are 



evaluated on the background (f^^° = and = Q,„). 

We evaluated the magnitude of the fluctuations in Dy, based 
on the cosmological parameters of the concordance model, as 
presented in Fig. 9. Fluctuations as low as one per cent are 
reached for much smaller domains than for the cosmic variance 
of the Q-parameters. Thus, at first sight it might seem that the 
BAO measurement of Dy could essentially overcome the cosmic 
variance limit. Closer inspection of this result reveals that this is 
not the case. Indeed, the much smaller variation of the distance 
Dy means that a precise knowledge of the distance measure Dy 
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Figure 9. Errors on the distance Dy for various survey geome- 
tries as a function of maximum redshift. For comparison the er- 
ror induced by the finite number of BAO modes in the corre- 
sponding full sphere volume, calculated with the fitting formula 
of Seo & Eisenstein (2007), is shown (insufficient volume). This 
error is about a factor of 10 bigger than the error from the local 
volume distortion caused by inhomogeneities that we calculated. 
Adding a shot noise term, corresponding to a galaxy density of 
« = 3 X 10 "^/z^Mpc^^ typical for SDSS and BOSS, we find that 
the cosmic variance of Dy is a subdominant contribution to the 
error budget. 



does not lead to an equivalently good estimate of the cosmic pa- 
rameters. 

Clearly, the systematic uncertainty that we calculated is only 
a minor effect compared with the eiTors intrinsic to the actual 
measurement of the acoustic scale, as a comparison of the three 
solid (green) lines in Fig. 9 shows. The lowest one is the fluctua- 
tion of the scale Dy for full spheres of the corresponding size at 
different places in the universe. It is therefore the possible local 
deformation caused by statistical over- or underdensities. The 
possible precision of a measurement of Dy by BAOs, however, 
also depends on the number of observable modes. This induces 
an error if the volume is too small, and in particular when it 
is smaller than the BAO scale a reasonable measurement is no 
longer possible. Accordingly, even for a perfect sampling of the 
observed volume, the error will not be smaller than the solid 
(green) lines in the middle. If one adds shot noise caused by im- 
perfect sampling by a galaxy density of « = 3 x 10^'^/i^Mpc"^, 
typical for SDSS and BOSS, the error increases even more. This 
means that for the realistic situation where we do not have a suf- 
ficiently small perfect ruler to allow for large statistics akeady 
for the small volumes considered here, the deformation uncer- 
tainty that we calculated remains completely subdominant. 

8. Fluctuations of the Hubble scale 

Local fluctuations of the Hubble expansion rate have already 
been considered in the literature (Turner et al., 1992; Shi et al., 
1996; Wang et al., 1998; Umeh et al., 2010). Here we wish to 
add two new aspects. 

The first one is on the measurement of H{z) itself. 
Experiments that try to measure i/ as a function of z, like the 
WiggleZ survey (Blake et al., 201 1), do this by measuring a "lo- 
cal" average H(z,n) in a region around the redshift z,„. These re- 
gions should not be too small to keep the effects of local fluc- 
tuations small. On the other hand they cannot be enlarged in an 



Figure 10. Optimum thickness of shells to minimize the variance 
of H{z) (see text for the two competing effects). The respective 
error corresponds to the height of the bars, the shells necessary 
for this purpose to their width. 



arbitrary way because then the redshift z„, becomes less and less 
characteristic for the averaging domain. In other words, for an in- 
creasingly thicker shell Az, the evolution of Hiz) begins to play 
a role. Therefore, one may find the optimal thickness of the av- 
eraging shells over which the variation in the expansion rate 



Var[i/(z)] 



Vvl 



H[zir)fWo(r)d'r 



H[{z(r)]Wo{r)d'r 



(65) 



equals the variance imposed by the inhomogeneous matter distri- 
bution. The corresponding shells are shown in Fig. 10. It should 
be noted that the error for the first bin is certainly underesti- 
mated in our treatment, which rests on linear perturbation the- 
ory. Taking into account higher orders, which become dominant 
at small scales, will certainly increase it. Of course, in these mea- 
surements the survey geometries will not necessarily be close to 
the SDSS or the 2dF geometry, but they are shown to illustrate 
survey geometries that do not cover the full sky. 

Secondly, we wish to note that the relation between fluctua- 
tions in the Hubble expansion rate and fluctuations in the matter 
density offers the interesting possibility to determine the evo- 
lution of the growth function for matter perturbations from the 
variances of the Hubble rate measured at different redshifts. A 
direct measurement of the growth function by a determination 
of erg at different epochs is difficult, because one never exam- 
ines the underlying dark matter distribution. Therefore one has 
to assume that the observed objects represent the same cluster- 
ing pattern as the underlying dark matter (this is the problem of 
bias). It is well known that there is bias and its modeling typi- 
cally has to rely on assumptions. 

An interesting bypass is to look at the variation of local ex- 
pansion rates at different redshifts. The assumption that the lu- 
minous objects follow the local flow is more likely and the as- 
sumption that this local flow is generated by the inhomogeneities 
of the underlying dark matter distribution is also reasonable. A 
similar idea leads to the attempt to use redshift-space distortions 
to do so (Percival & White, 2009). The fact that one considers 
fluctuations means that we would not have to know the actual 
value of H(z), but only the local variation at different redshifts. 
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This variation, defined as 

Ho - H^fiao) 
Oh - — = , 

has the fluctuations of Eq. (47) 

o-{_5h) = ^^^2)(ao)/£)(ao)cro„ . 



(66) 



(67) 



If we were to measure this quantity at different redshifts, we 
could, without knowledge of the absolute normalization of 
//© (z), determine (a®) only from the variance and therefore 
the constant c - Q./^/Q.,„. 

Note that in the standard case, where the background redshift 
is identified with the observed one, (a©) is simply replaced 
by the growth rate f(a) - and measuring the Hubble 

fluctuations would yield a direct measurement of /. In the real 
world, where the redshift captures the structure on the way from 
the source to us, it is not dkectly the background redshift. One 
would rather measure the modified "growth rate" (flo). The 
difference between these two quantities is small in our range of 
validity for (a©), however (corrections of linear order in the 
perturbations). 



9. Conclusion 

For the first time, we brought together the well-established per- 
turbative approach to incorporate inhomogeneities in Friedman- 
Lemaitre models (the theory of cosmological perturbations) and 
the ideas of cosmological backreaction and cosmic averaging in 
the Buchert formalism. Focusing on observations of the large- 
scale structure of the Universe at late times, we showed that the 
cosmic variance of cosmological parameters is in fact the lead- 
ing order contribution of cosmological averaging. 

We studied volume averages, their expected means, and vari- 
ances of the cosmological parameters Hq, Q^, Q^, Qa (Q^ is of 
higher order in perturbation theory). Our central result is sum- 
marized in (47)-(51). 

Our extension of the backreaction study of Li & Schwarz 
(2008) to the ACDM case enabled us to study fluctuations for 
a wider class of cosmological models. We were able to confirm 
for the fluctuations in the Hubble rate that our results in comov- 
ing synchronous gauge agree with those found in Poisson gauge 
(Umeh et al., 2010). 

The use of general window functions allowed us to consider 
realistic survey geometries in detail and to calculate the fluctua- 
tions in the matter density, empirically found in the SDSS data 
by Driver & Robotham (2010), directly from the underlying DM 
power spectrum. Converting this information into curvature fluc- 
tuations, we found that regions of 540/;" 'Mpc diameter may still 
have a curvature parameter of ~ 0.01, even if the background 
curvature vanishes exactly. We found that cosmic variance is a 
limiting factor even for surveys of the size of the SDSS survey. 
A volume-limited sample up to a redshift of 0.5 was able to con- 
strain the local curvature to 0. 1 per cent. 

Finally, we investigated the distortions of the local distance 
to a given redshift and found that it is less affected by the fluc- 
tuations of the local cosmic parameters than one might expect. 
The distance measure Dy, used in BAO studies, is accurate to 
0.2 percent for samples ranging up to z a; 0.35. This means that 
BAO studies are not limited by cosmic variance, but by problems 
such as sampling variance and insufficient volume, as discussed 
in section 7. 



In a next step one should incorporate the second-order ef- 
fects into the expected means of the cosmological parameters. 
There are no second-order corrections to the variances, as ar- 
gued in section 2. Therefore a complete second-order treatment 
seems feasible. 

The limitation of our approach comes from the fact that 
Buchert's formalism relies on spatial averaging. Averaging on 
the light cone would be more appropriate (Gasperini et al., 
20 11), thus the study in this work has been restricted to redshifts 
<K 1, where we expect light cone effects to play a subdominant 
role. 
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